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Abstract. The Sunyaev-Zeldovich (fs^; SZ) effect has proven to be an extremely powerful 
tool to study the physical and evolutionary properties of rich clusters of galaxies. Upcoming 
SZ experiments, with their much improved sensitivity and angular resolution, will provide 
unique information also on phases of galaxy evolution characterized by the presence of 
large amounts of hot proto-galactic gas. We present a preliminary analysis of the SZ signals 
that can be expected at the collapse of the proto-galaxy when, according to the standard 
scenario, the gas is heated at its virial temperature, and during episodes of strong energy 
injections from the active nucleus. The contributions of such signals to excess power on 
arc-minute scales recently found by CBI and BIMA experiments are briefly discussed. 

INTRODUCTION 

The currently standard hierarchical clustering paradigm for large scale structure 
formation in a ACDM universe has successfully confronted a broad variety of 
observations, ranging from the formation of galaxy clusters, to large scale velocity 
fields, to power spectra of the galaxy distribution and of the microwave background, 
and more (see, e.g., jH^). However, the theory of formation and evolution of galaxies 
is not in a very satisfactory state. Serious challenges have emerged in the last 
years (jH^]): the excess of predicted small scale structure; the persistent inability 
of even the best semi-analytic models ( [HI [13, ) to account for the surface 

density of massive galaxies at substantial redshifts detected by (sub)-mm surveys 
with SCUBA and MAMBO (|1[I3) and by deep K-band surveys ([la^]); the 
difficulties to account for the distribution of velocity dispersions of low-ionization 
damped Lyman-a systems at ^ > 1.5 ([il,[i^); the low predicted specific angular 
momentum of galactic disks (0,[i3|); the observational evidence contradicting the 
existence of the central cusp in the dark matter distribution, predicted by numerical 
simulations ([13, [I^ [i^). On larger scales, the observed relationship between X- 
ray luminosity and gas temperature for groups of galaxies strongly deviates from 
expectations of the simplest hierarchical clustering models. 

Solutions of these problems may require advances in different fields that may 



include particle physics (self-interacting dark matter?), deviations from a power law 
of the power spectrum of primordial density perturbations (as may be suggested 
by WMAP data, ^^J), or the astrophysics of galaxy formation and evolution. In 
any case, a better understanding of the complex physical processes governing the 
galaxy formation is mandatory. 

We discuss here how unique information may be provided by the Sunyaev- 
Zeldovich (SZ) effect. In fact, the the proto-galactic gas is expected to have a 
large thermal energy content, leading to a detectable SZ signal, both when the 
protogalaxy collapses with the gas shock- heated to the virial temperature (j50,EH), 
and in a later phase as the result of strong feedback from a flaring active nucleus 

(|ilHHIiil,.42,izj). 



GALAXY-SCALE SUNYAEV-ZELDOVICH EFFECT 

Let us consider a fully ionized gas with a thermal energy density egas, within the 
virial radius 

where h is the Hubble constant in units of lOOkms"^ Mpc~^ and Pvir — 200p„, p„ = 
1.88 10^^^/i^(l + 2vir)^gcm~^ is the mean density of the universe at the virialization 
redshift ^vir- 

The Comptonization parameter y, characterizing the amplitude of the Sunyaev- 
Zeldovich effect due to this gas, can be estimated as ([g^): 

y^j- — ; (2) 

here ecMB = «bs^cmb — 4-2 x 10~^^(l + z)'^ergcm~^, gbb being the black-body con- 
stant and TcMB = Tq{1 + z) the temperature of the cosmic microwave background 
(CMB). The amount Ae of gas thermal energy transferred to the CMB through 
Thomson scattering is : 

(3) 

tc being the gas cooling time by Thomson scattering. 

A useful reference value for the thermal energy content of the gas, i^gas is its 
binding energy (Eb,gas = Mgas^^^i^, ^^vir = I62h'/^{1 + zy/^{M,iJ10'^M^YKms-' 
being the circular velocity of the galaxy at its virial radius, cf. |43|). The 
amplitude of the SZ dip in the Ray leigh- Jeans region can be written as: 

|AT|,, = 2,TcMB - l.T ( AWl±^)^^.a./M.. M 

' '^^ ^ \0.5j \ 3.5 J 0.1 lOi2M0Eb,gas 

Although current estimates of the ratio of the mass in stars to the total mass in 
massive spheroidal galaxies yield average values of Mgtar/^vir — 0.03 



FIGURE 1. Counts of galaxy-scale SZ effects due to proto-galactic gas heated at the virial 
temperature (dashed curve) and to quasar-driven blast- waves (solid), at 30 GHz. The flux S is 
the absolute value of ^sz [eq. ifTUIl ] 

value of Mgas/Mvir — 0.1, close to the cosmic ratio between baryon and total mass 
density, is likely to be more appropriate for proto-galaxies before a large fraction 
of the initial gas mass is swept off by the combined action of supernova explosions 
and of quasar feedback. 

The SZ effect in Eq. Q shows up on small (typically sub-arcmin) angular scales. 
Quasar-driven blast-waves could inject into the ISM an amount of energy several 
times higher than the gas binding energy, thus producing larger, if much rarer, SZ 
signals. 

SZ EFFECT FROM QUASAR FEEDBACK 

Impressively many lines of evidence converge in indicating that the formation of 
super-massive black holes powering nuclear activity is intimately linked to the 
formation of their host galaxies ([l7L l6^. These include: 

• the discovery that Massive Dark Objects (MDOs), with masses in the range 
~ 10^-3 X 10^ Mq and a mass function matching that of baryons accreted onto 
black holes during the quasar activity ([EH), are present in essentially all local 
galaxies with a substantial spheroidal component (pollscl. IHoj]): 

• the tight correlation between the MDO mass (Mr h) and the velocity dispersion 
(a) of stars in the host galaxy ([3, Is^. l58| ) . the mass of the spheroidal 
component ([3, E^), and the mass of the dark halo [13]; recently j53|] have 
found that the Mbh^c" correlation is already present at redshift up to ~ 3; 

• the correspondence between the luminosity function of active star-forming 
galaxies at 2; ~ 3, the B-band luminosity function of quasars, and the mass 
function of dark halos at the same redshift (j^H^I); 



• the similarity between the evolutionary histories of the luminosity densities of 
galaxies and quasars (e.g. j^). 

As a consequence, the evolutionary histories of both active nuclei and spheroidal 
galaxies (or galactic bulges) can onlv be understood if we properly allow for their 
mutual feedback. As discussed by ji^,!!^!), the mutual feedback between galaxies 
and quasars during their early evolutionary stages may indeed be the key to 
overcome some of the crises of the currently standard scenario for galaxy evolution. 

One ingredient of the scenario proposed by |23] are powerful quasar driven 
outflows, carrying a considerable fraction of the quasar bolometric luminosity, Lboi, 

1 /2 

increasing as Lj^^j (for emission close to the Eddington limit). For large enough 
galaxies (Mvir>lO^^M0), the interstellar medium (ISM) is eventually swept out, 
and, correspondingly, the star-formation is switched off, on shorter timescales for 
more massive objects. The outflows are expected to be highly supersonic, and 
therefore liable to induce strong shocks, transiently heating the interstellar gas to 
temperatures exceeding the virial value. 

A black-hole (BH) accreting a mass Mbh with a mass to radiation conversion 
efficiency esn releases an energy Ebu = eeH^BHC^. We adopt the standard value 
for the efficiency e^n = 0.1 and assume that a fraction fh = 0.1 of the energy is 
fed in kinetic form and generates strong shocks turning it into heat. Actually the 
fraction, fh, of energy going into heating of the interstellar gas may be a function 
of the BH mass, or of the bolometric luminosity of the active nucleus (AGN). For 
example, [l^l argue that, based on the model for AGN-driven outflows by jio], for 
emission at the Eddington limit, the fraction of bolometric luminosity released in 
kinetic form increases as lI/J. This conclusion is, however, model-dependent and 
we prefer the simpler assumption of a constant value for fh, that may be viewed 
as an effective, luminosity weighted value. As discussed by 0, 113], ~ 0.1 can 
account for the pre-heating of the in terg alactic medium in groups of galaxies. 

Using the recent re-assessment by jSSf of the well known correlation between the 
BH mass and the stellar velocity dispersion 

Mbh = 1.410^ ( ^—r-) Mq . (5) 

V200km/sy ° ^ ^ 



we get 
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(6) 



The amplitude of the SZ dip in the Rayleigh- Jeans region due to quasar heating 
of the gas is then estimated as: 
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FIGURE 2. Power spectra of galaxy-scale SZ effects at 30 GHz, compared to the CMB power 
spectrum and to the CBI (js^l; large square) and BIMA (E^! Point with 68% confidence error 
bars and 95% confidence upper limit) data. The dashed lines correspond to SZ effects due to proto- 
galactic gas, the solid lines to quasar-driven blast-waves. For each set, the steep straight line refer 
to Poisson fluctuations, the flatter line to fluctuations due to clustering, the uppermost line to 
the sum (in quadrature) of the two contributions. The dotted line is the sum (in quadrature) of 
all contributions. The downward arrows signify that, as mentioned in the text, the contributions 
of proto-galactic gas (dashed lines) are actually upper limits. We also argue that the CBI results 
(large square) might be, conservatively, interpreted as upper limits 



Following |47[, we adopt an isothermal density profile of the galaxy. The virial 
radius, encompassing a mean density of 200 pu, is then: 



corresponding to an angular radius: 



where dA^z) is the angular diameter distance. 



(8) 



/Ebh \ f eBH \ fl + z^\ ^^(2.5) 



COUNTS AND SMALL-SCALE FLUCTUATIONS 

The "flux density", which in the Rayleigh- Jeans region is actually negative, asso- 
ciated to the SZ effect is given by (neglecting relativistic corrections): 

2(A;Tcmb)^ / n /.^x 
Ssz{v) = — — 9{x)yuj , (10) 



where oj = tiO^z and 



with X = hu/kTcMB- 



x^e^ / e^' + 1 
T? V e^-1 



9{^)= ... .^2 (^-^-4) , (11) 



Quasar-driven SZ signals 

The counts of quasar-driven SZ signals, and the corresponding small scale 
fluctuations, can be estimated from the evolving B-band luminosity function of 
quasars, (f){LB,z), relating the total amount of energy released, E^n, to the B- 
band luminosity. Lb- If the quasar radiates at the Eddington limit E^u = ^boi-^B^Si 
where fcboi is the bolometric correction and ts — 4.5 x 10^ yr is the quasar lifetime 
for eBH = 0.1. 

The lifetime, tsz, of the transient SZ effect due to the quasar-driven blast-waves 
is approximately the time for the shock to reach the outer border of the host galaxy. 
Using the expression (j^) for the evolution with time of the radius of a self-similar 
blast-wave expanding in a medium with an isothermal density profile, p oc r~^, we 
have: 



0.5 / V 1062erg/ VO.l/ V 0.1 / V 3.5 



The number density of SZ sources per unit interval of SZ "fiux" S^z can then be 
estimated as ^ 

(psziSsz, z) = (f){LB,z)-^—^ , (13) 

where Lb{Ssz,z) is the blue luminosity of a quasar at redshift z causing an SZ fiux 
5'sz, and tq,opt is the duration of the optically bright phase of the quasar evolution. 

Figure 1 shows, as an example, the predicted counts at 30 GHz, as a function of 
5'sz, obtained adopting the exponential model for the evolving luminosity function 
of quasars derived by |4^ for an Einstein-de Sitter universe, an Hubble constant of 
50kms~^ Mpc~^, and an optical spectral index of quasars a = 0.5 {S„ oc z/"). The 
parameters have been set at een = 0.1, = 0.1, fcboi = 10, tq,opt = 10^ yr. 



SZ signals from proto-galactic gas at virial temperature 

As discussed by jl^ quasars can be used as effective signposts for massive 
spheroidal galaxies in their early evolutionary phases; this allows us to bypass 
the uncertainties in the normalization of the primordial perturbation spectrum af- 
fecting estimates based on the various versions of the Press & Schechter formalism. 
According to [1§], the mass of their dark-matter halo, Mvir, is related to the mass 



of the central black-hole by: 



Mbh M^ir ' 



0-1 T7^ • (14) 



Under the assumption that quasars emit at the Eddington limit, the number density 
of sources with gas at virial temperature, producing SZ dips of amplitude given by 
Eq. (jU, can be straightforwardly related to the quasar luminosity function (f){LB,z) 
by an equation analog to Eq. (|T!?jl . replacing tsz with the gas cooling time, tcooi- 

The counts shown in Fig. 1 refer to the same cosmological model and the same 
lifetime of the quasar optically bright phase as above and to Mgas/Mvir = 0.1. As for 
the cooling time, which is rather uncertain, we have made the extreme assumption 
that tcooi = tcxp- Thus the counts plotted have to be taken as upper limits. 



Small scale fluctuations 

The Poisson fluctuations are straightforwardly computed from the counts; since 
the latter are very steep, the results are insensitive to the upper flux density cutoff. 

To estimate fluctuations due to clustering we adopted the spatial correlation 
function of quasars estimated by for an Einstein-de Sitter universe (^(r) = 
(r/ro)~^'^^, with tq = 4.29/i~^ Mpc, constant in comoving coordinates), cut-off at 
r = Svq. 

The power spectra of fluctuations due to both the gas in virial equilibrium 
and heated by the quasar feedback are compared, in Fig. 2, with the recent 
measurements of arcminute scale fluctuations at 30 GHz by the CBI ([s^]) and 
BIMA (^3|) experiments. We stress again that our estimated of signal from the 
gas at the virial temperature is based on an over-simplified treatment. In fact, 
the thermal history of the protogalactic gas is quite complex: only a fraction of it 
may be heated to the virial temperature (0,01); cooling may be relatively rapid 
in the densest regions; on the other hand, significant heating may be provided 
by supernovae and by the central active nucleus. Also, we expect that (jl^), 
particularly in the case of the most massive halos virializing at high z, the quasar 
feedback removes the residual interstellar gas on a timescale ranging to 0.5 to a 
few Gyrs. Thus, again, the results shown in Fig. 2 are (possibly generous) upper 
limits. 



DISCUSSION AND CONCLUSIONS 

Evidences of statistically significant detections at 30 GHz of arcminute scale 
fluctuations in excess of predictions for primordial anisotropies of the cosmic 
microwave background (CMB) and also in excess of the estimated contamination 
by extragalactic radio sources have recently been obtained by the CBI (^]) and 
BIMA (^]) experiments. 



It may be noted that, although the CBI group apphed a quite elaborated 
treatment of radio- sources, the sensitivity of their 31 GHz point source observations 
does not really guarantee proper allowance for this foreground. Their strategy 
comprised pointed 31 GHz observations with a 4(T detection limit of Ssighz = 6mJy, 
of all NVSS sources with 5'i.4ghz > 6mJy and a direct survey reaching a flux limit 
of 5'iim,3iGHz = 6mJy, over an area of 1.8 square degrees. If the slope of their 
31 GHz counts keeps constant down to flux densities several times lower than 
>5'iim,3iGHz, the Poisson fluctuations due to sources below the detection limit amount 
to ~ 30/iK, and may entirely account for the measured signal. Thus, their excess 
signal (15-30 /iK in the multipole range ^ = 2000-4000) may be more conservatively 
interpreted as an upper limit. 

The BIMA group (jl3]) carried out a VLA survey at 4.8 GHz of their fields down 
to a flux density limit of ~ 150/iJy to identify and remove point sources. In this 
case, the residual contamination is indeed likely to be small, at the cost of a loss 
of sensitivity. 

If indeed the detected signal cannot be attributed to extragalactic radio sources, 
its most likely source is the thermal SZ effect (J2l|). The SZ within rich clusters of 
galaxies has been extensively investigated (jHEJ)- The estimated power spectrum 
was found to be very sensitive to the normalization (cts) of the density perturbation 
spectrum. A normalization erg > 1, somewhat higher than indicated by other 
data sets, is apparently required to account for the CBI data. High resolution 
hydrodynamical simulations of structure formation have highlighted the 

presence of substantial SZ structure on sub-arcmin scales. 

We have carried out an analytical investigation of the SZ signals associated to 
the formation of spheroidal galaxies. We find that they are potentially able to 
account for the BIMA results. Proto-galactic gas heated at the virial temperature 
and with a cooling time comparable with the expansion time may provide the 
dominant contribution; in this case we expect clustering fiuctuations of amplitude 
comparable to Poisson fiuctuations, although with a fiatter power spectrum. 

Potentially detectable SZ signals can also be produced by strong feedback from 
the central active nucleus. It is, in fact, becoming increasingly clear that the 
interplay between star-formation and nuclear activity plays a key role in shaping 
the early evolution of both the host galaxy and the active nucleus. The enormous 
amounts of energy that the active nuclei (AGNs) release to their environment, not 
only in radiative but also in kinetic form, can strongly infiuence the gas both in the 
host galaxy and in the surrounding intergalactic medium, with conspicuous effects 
on X-ray emission and on the SZ signal. A flaring AGN can drive an energetic shock, 
causing a transient, strong rise of the average gas pressure and, consequently, of 
the amplitude of the SZ effect. These signals are much rarer but can individually 
reach substantially higher values. 

High frequency surveys with sub-arcmin resolution can directly detect individual 
SZ sources. We expect a surface density of ~ ldeg~^ at S'sqghz — ImJy. Of course, 
the counts at these relatively bright fluxes are likely dominated by SZ effects in 
clusters of galaxies, which however should be easily distinguishable because of their 
much larger angular size and much lower redshift. 

We thus conclude that the SZ effect is an effective probe of the thermal state 



of the gas, of its evolution and, in particular, of the AGN feedback. The main 
parameter governing the global effect on the ambient medium of the energy released 
by AGNs is the ratio, AE/E, between the amount of energy deposited in the 
medium and the binding energy of the receiving gas (^8]). lS.EjE can plausibly 
be > 1 for very massive galaxies and small groups, but is -C 1 for rich clusters. 
Qualitatively different effects are therefore expected in the two cases. If AE / E>1 
the energy injection heats the gas and partially ejects it out of the potential 
well. This may lead directly to a correlation between the black-hole mass and the 
velocity dispersion of stars (or the bulge mass) with slope and normalization in nice 
agreement with the data. Also, the ensuing increase of the gas temperature and 
decrease of the X-ray luminosity (because of the decreased gas density) may explain 
the deviation of the observed T-Lx from predictions of the hierarchical clustering 
scenario, in the presence of gravitational effects only. Ongoing and forthcoming high 
resolution surveys in the cm/mm region may directly test this scenario. Clearly the 
SZ effect is a much more effective probe of the thermal energy content of the plasma 
than the X-ray emission. The former is in fact an almost perfect calorimeter (0), 
being directly proportional to the thermal energy. On the contrary. X-ray emission 
is very sensitive to the details of the density distribution of the emitting gas. 
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